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Abstract 

Background: The genome-wide single-nucleotide polymorphisms (SNPs) profiles can be used as diagnostic markers for 
human cancers. The associations between mitogen-activated protein kinase kinase kinase 1 (MAP3K1) SNPs rs889312 A>C, 
rsl 68861 65 T>G and breast cancer risk have been widely evaluated, but the results were inconsistent. To derive a 
conclusive assessment of the associations, we performed a meta-analysis by combining data from all eligible case-control 
studies up to date. 

Methods: By searching PubMed, ISI web of knowledge, Embase and Cochrane databases, we identified all eligible studies 
published before September 2013. Odds ratios (ORs) with 95% confidence intervals (CIs) were used to assess the strength of 
associations in fixed-effect or random-effect model. False-positive report probability (FPRP) was calculated to confirm the 
significance of the results. 

Results: A total of 59670 cases in 20 case-control studies were included in this meta-analysis. Significant associations with 
breast cancer risk were observed for SNPs rs889312 and rs16886165 polymorphisms with a per-allele OR of 1.1 1 (95% CI: 
1.09-1.13) and 1.14 (95% CI: 1.09-1.20) respectively. For rs889312, in subgroup analysis by ethnicity, significant associations 
were identified in Europeans and Asians, but not in Africans. When stratified by estrogen receptor (ER) expression status, 
rs889312 was associated with both ER-positive and ER-negative breast cancers. Results from the FPRP analyses were 
consistent with and supportive to the above results. 

Conclusions: The present meta-analysis suggests that rs889312-C allele and rsl 68861 65-G allele might be risk factors for 
breast cancer, especially in Europeans and Asians. 
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Introduction 

Globally breast cancer is the most commonly diagnosed cancer 
and the leading cause of cancer death in females, with an 
estimation of 1 million new cases and over 400,000 deaths per year 
[1]. Breast cancer is a genetically heterogeneous disease following 
a polygenic mode of inheritance. A better understanding of the 
genomic variations pertaining to the disease will eventually lead to 
improved diagnostic and treatment strategies for breast cancer 
patients. Evidence has been accumulated that mutations in several 
high and moderate penetrance genes, such as BRCA1, BRCA2, 
ATM, CHEK2, BRIP1 and PALB2, contribute to breast cancer 
susceptibility. However, these mutations explain only a small 
portion of the disease risk and the majority of genetic variations 
that contribute to the breast cancer risk remains unclarified, 
especially those of low penetrance genes. Genome-Wide Associ- 



ation Study (GWAS), an examination of disease-related genetic 
variations across the human genome, provides a powerful tool in 
identifying such mutations/polymorphisms genome-wide. Some 
common low penetrance genetic variants that might confer 
increased risk to breast cancer have been identified recently [2,3]. 

Mitogen-activated protein kinase kinase kinase 1 (MAP3K1) gene, 
which is located in the chromosome 5qll.2, encodes a serine/ 
threonine kinase which is involved in the mitogen-activated protein 
kinase (MAPK) signaling pathway. MAPK signal transduction 
regulates the transcription of important cancer genes including c- 
Myc, c-Elkl, c-Jun and c-Fos [4,5]. The single-nucleotide polymor- 
phism (SNP) rs889312 in MAP3K1 was identified to be associated 
with breast cancer risk by GWAS [5], with confirmation of the 
association in European ancestry population by another study [6] . A 
GWAS conducted in subjects of European descendant reported that 
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Records excluded 
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Full-text articles 
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Full-text articles 
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3 meta-analyses 
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(meta-analysis) 
(n = 20) 



Figure 1. Flowchart of the study selection process. 

doi:10.1371/journal.pone.0090771.g001 

the minor allele of MAP3K1 rs 16886 165 was associated with 
increased risk of breast cancer under heterozygote codominant and 
homozygote codominant genetic models [7]. 

Up to date, a number of studies have investigated the relationship 
between polymorphisms in MAP3K1 and breast cancer or its 
malignant phenotypes, whereas the results were inconsistent and 
inconclusive in populations of diverse ethnicities, especially among 
Asians and Africans [8,9]. The current GWASs were mainly 
conducted in people of European ancestry, and most single studies 
with insufficient sample sizes were not powerful in detecting the 
minor effects of low penetrance alleles on breast cancer risk. 
Therefore, the meta-analysis presented here aims to summarize the 
available evidence for the genomic variants in NIAP3K1 and intends 
to provide the highest level of evidence for the association between 
breast cancer risk and MAP3K1 (rs889312 and rsl6886165) 
polymorphisms among diverse ancestry populations. 

Materials and Methods 

Literature search strategy 

We conducted a comprehensive literature search in Pubmed, 
ISI web of knowledge, Embase and Cochrane databases up to 
September 2013 using the following terms: "breast cancer" or 
"breast carcinoma" and "A4AP3KF' or "mitogen-activated protein 
kinase kinase kinase 1" or "rs889312" or "rsl6886165" and 
"polymorphism" or "variation". All potentially relevant studies 
were retrieved, and the abstracts were screened to exclude clearly 
irrelevant studies. The remaining articles were read to determine 
whether they contained information on the topic of interest. We 
also examined the reference lists from the main reports and 
reviews to identify additional relevant studies. 



Selection criteria 

Eligible studies had to meet the following criteria: (1) original 
papers containing independent data which have been published in 
peer-reviewed journals, (2) case-control studies which evaluated 
the association between breast cancer risk and AIAP3K1 (rs889312 
and rs 16886 165) polymorphisms, (3) had odds ratio (OR) with its 
95% confidence interval (CI) or genotype distribution information 
for estimating OR (95% CI). 

Data extraction 

Two investigators independently reviewed and extracted data 
with any discrepancies resolved by consensus. For each eligible 
study, the following data were extracted: first author's surname, 
year of publication, country of origin, population ethnicity, source 
of control subjects, genotyping method, age, the values of OR with 
its 95% CI, total numbers of cases and controls, the genotype 
counts in cases and controls, matching criteria, whether genotype 
distribution of control group was consistent with Hardy- Weinberg 
equilibrium (HWE). Ethnic groups were categorized as European, 
Asian and African (i.e. people of European, Asian and African 
ancestry). Different case-control groups in one study were 
considered as independent studies if the data was available. 

Quality assessment 

Methodological quality was assessed using a classification 
method known as "extended-quality score". The criteria cover 
diagnostic criteria, the degree of matching of controls, genotyping 
examination, Hardy-Weinberg equilibrium in the control popu- 
lation, and bias in data processing. Each paper was scored of 
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Study 
ID 



OR (95% CI) 



% 

Weight 



European 

Garcia-Closas(2008) 

Rebbeck(2009) 

Barnholtz-Sloan(2010) 

Gorodnova(2010) 

Latif(2010) 

Tamimi(2010) 

Campa(2011) 

Campa(2011) 

Campa(2011) 

Campa(2011) 

Campa(2011) 

Slattery(2011) 

Butt(2012) 

Harlid(2012) 

Shan(2012) 

Subtotal (l-squared = 0.0%, p = 



0.708) 



Asian 

Garcia-Closas(2008) 

Zheng(2010) 

Campa(2011) 

Han(2011) 

Jiang(2011) 

Chan(2012) 

Kim(2012) 

Subtotal (l-squared = 42.0%, p = 0.1 1 1 ) 
African 

Rebbeck(2009) 
Barnholtz-Sloan(2010) 
Campa(2011) 
Chen(2011) 

Subtotal (l-squared = 40.2%, p = 0.171) 
NA 

Jara(2013) 
Resler(2013) 

Subtotal (l-squared = 0.0%, p = 0.442) 

Heterogeneity between groups: p = 0.016 
Overall (l-squared = 27.0%, p = 0.095) 



T 




1.12(1.09, 1.15) 
1.10(0.90, 1.35) 
1.18(1.03, 1.34) 

1.01 (0.71, 1.44) 
1.17(0.96, 1.41) 
1.12(0.95, 1.31) 
1.21 (1.02, 1.43) 
1.13(1.04, 1.23) 
1.11 (1.02, 1.21) 
1.06 (0.89, 1.26) 
1.05 (0.93, 1.18) 

1.02 (0.93, 1.13) 
1.19(1.04, 1.37) 
1.13 (1.06, 1.21) 
1.33(1.09, 1.63) 
1.12(1.10, 1.14) 



1.12(0.97, 1.28) 
1.07 (0.99, 1.15) 
1.14(0.97, 1.36) 
1.22(1.06, 1.41) 
0.99 (0.83, 1.18) 
0.96 (0.85, 1.08) 
1.16(1.06, 1.27) 
1.09(1.05, 1.14) 



1.33 (0.98, 1.81) 
0.95 (0.80, 1.13) 
0.93 (0.75, 1.14) 
1.07 (0.99, 1.18) 
1.04 (0.97, 1.12) 



1.36(1.14, 1.63) 
1.24(1.06, 1.44) 
1.29(1.15, 1.45) 



1.11 (1.09, 1.13) 



43.70 

0.76 

1.81 

0.25 

0.85 

1.22 

1.10 

4.46 

4.30 

1.04 

2.21 

3.31 

1.65 

7.16 

0.77 

74.60 



1.63 
5.59 
1.10 
1.54 
1.01 
2.19 
3.84 
16.91 



0.33 
1.05 
0.72 
4.07 
6.17 



0.98 
1.34 
2.32 



100.00 



Figure 2. Forest plot of MAP3K1 rs889312 polymorphism and breast cancer risk stratified by ethnicity. Fixed-effect model was used for 
the analysis (allele contrast model C vs. A). 
doi:1 0.1 371 /journal.pone.0090771 .g002 



'high', 'median' or 'poor' quality. Detailed procedure of the 
quality assessment was performed according to literature [10]. 

Statistical methods 

The meta-analysis evaluated the association between rs889312 
and the risk of breast cancer, for the: (1) allele contrast model (C 



vs. A), (2) the heterozygote codominant model (AC vs. AA), and (3) 
the homozygote codominant model (CC vs. AA) model [1 1] . For 
rs 16886 165 polymorphism, pooled OR was obtained by an allele 
contrast model (G vs. T). The f-ttsX was performed to assess 
possible heterogeneity among studies [12]. Hardy- Weinberg 
equilibrium (HWE) was tested by the x 2 test. Pooled odds ratios 
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Study 



OR (95% CI) 



Weight 



Chen(2011) 



Chan(2012) 



Kim(2012) 



Zheng(2013) 



Overall (l-squared = 0.0%, p = 0.985) 



<i> 



1.15 (1.06, 1.25) 



1.15(1.01, 1.29) 



1.14 (1.04, 1.25) 



1.12 (1.00. 1.26) 



1.14 (1.09. 1.20) 



36.15 



16.41 



18.40 



100.00 



Figure 3. Forest plot of MAP3K1 rsl 68861 65 polymorphism and breast cancer risk. Fixed-effect model was used for the analysis (allele 

contrast model G vs. T). 

doi:1 0.1 371 /journal.pone.0090771 .g003 



(ORs) with 95% confidence intervals (95% CIs) were calculated 
using fixed-effects (Mantel-Haenszel method) [13] or random- 
effects (DerSimonian-Laird method) [14] models. When the effects 
were assumed to be homogeneous, the fixed-effects model was 
then used; otherwise, the random-effects model was more 
appropriate. Subgroup analyses were performed by ethnicity, 
sample size (<1000 cases and > 1000 cases) and estrogen receptor 
(ER) status (ER-positive and ER-negative). Sensitivity analyses 
were performed to evaluate the stability of the results by sequential 
omission of individual studies in the meta-analysis to show the 
influence of the individual data set to the pooled OR. The Egger's 
test and funnel plots were utilized to provide diagnosis of 
publication bias [15]. All of the above analyses were performed 
using ST ATA version 12.0 (StataCorp, College Station, TX). 

For each statistically significant association, the false-positive 
report probability (FPRP) analysis was performed using the method 
reported by Wacholder et al. [16]. A prior probability of 0.001 was 
set to detect an OR of 1.20. An FPRP cutoff value of 0.2 was used, 
and only the results with FPRP values less than 0.2 were referred as 
noteworthy. The Excel spreadsheet provided by Wacholder et al. 
was used to calculate statistical power and FPRP values (http:/ /jnci. 
oxfordjournals.org/ content/ 96/6/ 434/ suppl/DC 1). 

Results 

Study Characteristics 

Based on the above criteria, a total of 20 eligible studies 
involving 59,670 cases and 66,862 controls met all selection 
criteria for the pooled analyses [6,8,9,17-33], Figure 1 illustrates 
the study selection process and Table 1 describes the main 



characteristics of these studies. For the rs889312 polymorphism, 
18 studies including a total of 54674 cases and 64542 controls were 
included in the meta-analysis. For the rsl 6886 165 polymorphism, 
4 studies involved a total of 7949 cases and 7660 controls. Of the 
cases, 67% were Europeans, 21% were Asians, and 10% were of 
African descendants. The extended-quality scores ranged from 5 
to 8, and 19 were given high quality, while the remaining 1 study 
was given median quality. The distributions of genotypes in the 
controls in 18 studies were consistent with the Hardy- Weinberg 
equilibrium, but those in 2 other studies remained unknown. 

Meta-analysis results 

Overall, the rs889312 polymorphism was significantly associat- 
ed with increased breast cancer risk under all genetic models (C vs. 
A: OR= 1.11, 95% CI: 1.09-1.13, P=0.000, I 2 = 27.0%; AC vs. 
AA: OR= 1.09, 95% CI: 1.05-1.13, P= 0.000, 1 2 = 20.0%; CC vs. 
AA: OR=1.21, 95% CI: 1.15-1.28, P= 0.000, I 2 = 31.9%; 
Figure 2). Additionally, the meta-analysis resulted in a statistically 
significant association between rsl 6886 165 and breast cancer. The 
pooled OR for risk G allele was 1.14 under allele contrast model 
(95% CI: 1.09-1.20, P= 0.000, I 2 = 0.0%; Figure 3). 

In the subgroup analysis by ethnicity, our results indicated a 
significant association between the SNP rs889312 and breast 
cancer incidence in European ancestry populations (C vs. A: 
OR= 1.12, 95% CI: 1.10-1.14, P=0.000, I 2 = 0.0%; AC vs. AA: 
OR= 1.12, 95% CI: 1.08-1.18, P= 0.000, I 2 = 0.0%; CC vs. AA: 
OR= 1.24, 95% CI: 1.15-1.33, P= 0.000, I 2 = 0.0%). For Asians, 
two genetic comparisons produced significandy increased risks (C 
vs. A: OR= 1.09, 95% CI: 1.05-1.14, P= 0.000, I 2 = 42.0%; CC 
vs. AA: OR=1.15, 95% CI: 1.05-1.26, P=0.002, I 2 = 17.8%), 
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Table 2. Stratified analyses for MAP3K1 rs889312 polymorphism and breast cancer risk. 





Analyses 


Allele contrast model 






Heterozygote codominant model 




Homozygote codominant model 






OR (95% CI) 


P(Z) 


I 2 


Power 

OR, 1.2 FPRP 


OR (95% CI) 


P(Z) 


I 2 


Power 
OR,1.2 FPRP 


OR (95% CI) 


P(Z) 


I 2 


Power 
OR, 1.2 


FPRP 


Overall 


1.11(1.09-1.13) 


0.000 


27.0% 


1.000 


0.000 


1.09(1.05-1.13) 


0.000 


20.0% 


1.000 


0.003 


1.21(1.15-1.28) 


0.000 


31.9% 


0.386 


0.000 


Ethnicity 


European 


1.12(1.10-1.14) 


0.000 


0.0% 


1.000 


0.000 


1.12(1.08-1.18) 


0.000 


0.0% 


0.995 


0.020 


1.24(1.15-1.33) 


0.000 


0.0% 


0.180 


0.000 


Asian 


1.09(1.05-1.14) 


0.000 


42.0% 


1.000 


0.142 


1.02(0.95-1.11) 


0.560 


37.2% 


1.000 


0.998 


1.15(1.05-1.26) 


0.002 


17.8% 


0.819 


0.768 


African 


1.04(0.97-1.12) 


0.239 


40.2% 


1.000 


0.997 


0.95(0.80-1.12) 


0.544 


39.3% 


0.941 


0.998 


1.03(0.78-1.36) 


0.820 


51.3% 


0.859 


0.999 


NA 


1.29(1.15-1.45) 


0.000 


0.0% 


0.113 


1.148 


1.14(0.96-1.34) 


0.134 


0.0% 


0.733 


0.993 


1.86(1.42-2.43) 


0.000 


0.0% 


0.001 


0.891 


Sample size 


<1000 


1.12(1.07-1.16) 


0.000 


35.2% 


1.000 


0.000 


1.09(1.02-1.17) 


0.011 


27.7% 


0.996 


0.945 


1.32(1.19-1.47) 


0.000 


27.5% 


0.041 


0.010 


>1000 


1.11(1.09-1.14) 


0.000 


18.8% 


1.000 


0.000 


1.09(1.05-1.14) 


0.000 


14.7% 


1.000 


0.142 


1.18(1.11-1.26) 


0.000 


29.3% 


0.692 


0.001 



doi:1 0.1 371 /journal.pone.0090771 .t002 



but the heterozygote codominant model did not reach statistically 
significance. Among Africans, we did not detect any significant 
association under all genetic models (Table 2). According to 
sample size, the rs889312 polymorphism presented significantly 
increased risks of breast cancer both in small and large studies. 
The data on rs889312 polymorphism stratified by ER status were 
available in 6 studies involving 29,200 cases and 44,104 controls. 
Subsidiary analyses of ER status yielded a per-allele OR for 



ER-positive tumors of 1.12 (95% CI: 1.10-1.15, P=0.000, 
I 2 = 0.0%) and ER-negative tumors of 1.08 (95% CI: 1.04-1.12, 
P= 0.000, 1 2 = 29.6%; Figure 4). FPRP values at the pre-specified 
prior probability of 0.001 were lesser than 0.2, except associations 
among Asians under the homozygote codominant model and 
small group under the heterozygote codominant model (Table 2). 



Study 
ID 



OR (95% CI) 



% 

Weight 



ER+ 

Garcia-Closas 

Garcia-Closas 

Rebbeck 

Rebbeck 

Cam pa 

Chen 

Kim 

Shan 

Subtotal (l-squared = 0.0%, p = 0.718) 



ER- 

Garcia-Closas 

Garcia-Closas 

Rebbeck 

Rebbeck 

Cam pa 

Chen 

Kim 

Shan 



Subtotal (l-squared = 29.6%, p = 0.192) 

Heterogeneity between groups: p = 0.083 
Overall (l-squared = 14.2%, p = 0.291) 



o 

I 

0 



1.12 (1.09, 
1.10 (0.93, 

1.21 (0.96, 
1.57 (1.08, 

1.13 (1.07, 
1.09 (0.98, 
1.15 (1.04, 

1.22 (0.93, 
1.12 (1.10, 



1.08 (1.02, 
0.96 (0.76, 
0.69 (0.44, 
1.13 (0.66, 
1.03 (0.93, 
1.06 (0.94, 
1.18 (1.05, 
1.31 (0.99, 
1.08 (1.04, 



1.16) 45.37 
1.30) 1.57 
1.52) 0.83 
2.29) 0.31 
1.19) 15.56 
1.19) 4.66 
1.27) 4.40 
1.60) 0.60 
1.15) 73.29 



1.14) 14.21 
1.22) 0.78 
1.08) 0.22 
1.93) 0.15 
1.14) 4.24 
1.19) 3.16 
1.32) 3.36 
1.71) 0.59 
1.12) 26.71 



1.11 (1.09, 1.14) 100.00 



Figure 4. Forest plot of MAP3K1 rs889312 polymorphism and breast cancer risk stratified by ER status. Fixed-effect model was used for 
the analysis (allele contrast model C vs. A). 
doi:1 0.1 371 /journal.pone.0090771. g004 
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Begg's funnel plot with pseudo 95% confidence limits 

2 




-2 

H l 
0 .2 .4 .6 

s.e. of: logor 



Figure 5. Begg's funnel plot of MAP3K1 rs889312 polymorphism and breast cancer risk. 

doi:10.1371/journal.pone.0090771.g005 



Sensitivity analyses and publication bias 

A sensitivity analysis was performed by sequential removal of 
individual studies. The results suggested that no individual study 
significantly altered the pooled ORs. The publication bias of the 
studies was evaluated by funnel plot and Egger's test. As shown in 
Figure 5 and 6, the shape of the funnel plots did not reveal obvious 
evidence of asymmetry for both polymorphisms. All the p values of 
Egger's tests were more than 0.05, suggesting that no publication 
bias was detected in these studies. 



Discussion 

The common variants rs889312 and rsl6886165 lie in a linkage 
disequilibrium (LD) block of approximately 280 kb which includes 
MAP3K1 gene [31]. The MAP3K1 gene, which is also known as 
MEKK1 (Mek kinase 1), encodes a 196-kDa serine/threonine 
protein kinase that activates the ERK (extracellular signal- 
regulated kinase), JNK (c-Jun NH2-terminal kinase) and NF-kB 
(nuclear factor-KB) pathways [34]. The downstream signal 
transductions regulate the survival, differentiation, proliferation 
and apoptosis of cell, and appear to be involved in tumor 
development and tumor progression [35-37]. As rsl6886165 and 



Begg's funnel plot with pseudo 95% confidence limits 




s.e. of: logoi 



Figure 6. Begg's funnel plot of MAP3K1 rsl 68861 65 polymorphism and breast cancer risk. 

doi:10.1371/journal.pone.0090771.g006 
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rs889312 are 88.3 and 79.5 kb upstream of MAP3K1 respectively, 
the causal variant may be closer to rs 16886 165 than rs889312 [9]. 
Although these two SNPs do not change the structure and 
therefore the biochemical function of MAP3K 1 , they may have 
effects on the modulation of MAP3K 1 expression and therefore 
the tuning of MAPK signal transductions. 

This meta-analysis demonstrated that the rs889312 and 
rs 16886 165 SNPs in MAP3K1 were associated with increased 
breast cancer susceptibility. There are a number of meta-analyses 
that have been performed to elucidate the associations between 
SNPs and breast cancer risk in recent years; however, only 21.7% 
of those significant associations were noteworthy [38]. Therefore, 
we performed FPRP analyses in this work to confirm robustness of 
the linkages. Although a previous meta-analysis directed by Lu 
et al. [1 1] has demonstrated that the rs889312 polymorphism was 
significandy correlated with breast cancer risk, it only included 
seven case-control studies without subgroup analysis. In the 
current meta-analysis, we have made much more powerful and 
detailed analysis to support our results: (1) more studies were 
included; (2) subgroup analysis was conducted and stratified by 
ethnicity and sample size; (3) the association between ER status 
and breast cancer was considered. 

When stratified by ethnicity, the rs889312-G allele showed to be 
a risk factor for the development of breast cancer in European and 
Asian ancestry populations, but not in Africans. Different 
frequencies of mutant alleles in different ethnic groups may 
contribute to different susceptibilities to cancer [39]. Different life 
styles and environmental factors among the different populations 
lead to diverse gene-environment interactions, and may therefore 
also account for different cancer susceptibilities [40]. In addition, 
the result may be partly biased because of the limited number of 
studies in Africans, which had insufficient power to detect a 
different linkage disequilibrium (LD) pattern or minor effect of the 
SNP among African populations. Large population studies, which 
significandy reduce the evidence from smaller studies, play 
dominant roles in the meta-analysis. When considering sample 
size, we found rs889312 was associated with increased breast 
cancer risk both in large and small studies under all genetic 
models. For rs 16886 165, the available studies were insufficient to 
stratified analysis. The pooled result of the current meta-analysis 
indicated that the rsl6886165-G allele conferred breast cancer risk 
in Asian and African ancestry population. Additionally, 
rs 16886 165 have been identified as a low-penetrance risk factor 
for breast cancer in European ancestry population by GWAS [7] . 

The prognosis of breast cancer is affected by ER status. Previous 
studies suggested that rs889312 was related to ER-positive breast 
cancers, which was confirmed by our stratified analysis by ER 
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